Microcanonical Analyses of Peptide Aggregation Processes 
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We propose the use of microcanonical analyses for numerical studies of peptide aggregation transi- 
tions. Performing multicanonical Monte Carlo simulations of a simple hydrophobic-polar continuum 
model for interacting heteropolymers of finite length, we find that the microcanonical entropy be- 
haves convex in the transition region, leading to a negative microcanonical specific heat. As this 
effect is also seen in first-order-like transitions of other finite systems, our results provide clear evi- 
dence for recent hints that the characterisation of phase separation in first-order-like transitions of 
finite systems profits from this microcanonical view. 
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Thermodynamic phase transitions in macroscopic, in- 
finitely large systems are typically analysed in the ther- 
modynamic limit of a canonical ensemble, i.e., the tem- 
perature T is treated as an intensive external control pa- 
rameter adjusted by the heat bath, and the total sys- 
tem energy E is distributed according to the Boltzmann- 
Gibbs statistics. The probability for a macrostate with 
energy E is given by p{E) = g{E) exp{-E/kBT)/Z , 
where g{E) is the density of states, Z the partition sum, 
and ks the Boltzmann constant. As long as the mi- 
crocanonical entropy S{E) — ks^'ng^E) is a concave 
function of E, the microcanonical (caloric) temperature 
T{E) = [dS{E)/dE)-^ for fixed volume V and par- 
ticle number A/" is a monotonically increasing function 
of E. Consequently, the microcanonical specific heat 
Cv{E) = dE/dT{E) = -{dS/dE)^/{d^S/dE^) is posi- 
tive. The specific heat can only become negative in an 
energetic regime, where S{E) is convex. In this region, 
the caloric T{E) curve exhibits a typical backbending, 
which means that the system becomes colder with in- 
creasing total energy. For this reason, the temperature 
T is not the most appropriate control parameter and the 
analysis of such, in particular finite, systems is more ade- 
quately performed in the microcanonical ensemble, where 
the system energy E is considered as the adjustable ex- 
ternal parameter [l|, Q . 

It is a surprising fact that the backbending effect is in- 
deed observed in transitions with phase separation. Al- 
though this phenomenon has already been known for a 
long time from astrophysical systems Q, it has been 
widely ignored since then as somehow "exotic" effect. Re- 
cently, however, experimental evidence was found from 
melting studies of sodium clusters by photofragmenta- 
tion y. Bimodality and negative specific heats are also 
known from nuclei fragmentation experiments and mod- 
els @, @, as well as from spin models on finite lattices 
which experience first-order transitions in the thermody- 
namic limit 0, Q . This phenomenon is also observed in 
a large number of other isolated finite model systems for 
evaporation and melting effects d, ■ 

In this Letter, we demonstrate the usefulness of the 



microcanonical ensemble for studies of the aggregation 
process of small proteins (peptides), which, due to the 
fixed inhomogeneous sequence of amino acids, are neces- 
sarily systems of finite size. Understanding protein ag- 
gregation is essential not only for gaining insights into 
general mechanisms of protein folding, but also for un- 
raveling the reasons of human diseases caused by pro- 
tein clustering. A well-known example is associated with 
Alzheimer's disease, where a few identical small frag- 
ments of large proteins show the tendency to form fibrils, 
e.g., the hydrophobic A/?ig_22 segment of the /3-amyloid 
peptide A/3 

Our results are based on computer simulations of a 
simple continuum aggregation model for heteropolymers. 
Since the hydrophobic force governs the tertiary fold- 
ing process resulting in a compact hydrophobic core sur- 
rounded by a shell of mainly polar residues, in our model 
the 20 amino acids naturally occuring in proteins are 
classified as hydrophobic (A) and polar (B) [13]. For 
the individual peptides, we employ the AB model in 
three spatial dimensions. At a mesoscopic length scale, 
this coarse-grained model with virtual peptide bonds and 
virtual bond angles has proven quite successful in the 
qualitative characterisation of naturally observed protein 
folding channels [l4|. Keeping the same parameter sets 
for the interaction of monomers of different polymers, the 
model for the aggregate reads: 

^ = E 41] + E E *(''v.. ; '^v > ) - (1) 

where /i, v label the M polymers interacting with each 
other, and i^, index the N monomers of the /ith poly- 
mer whose intrinsic energy is given by 

4b = 4 E(l-^os^v)+ E 'J'(%J.;'^V''^jJ' 

(2) 

with < di^^ < IT denoting the bending angle between 
monomers i^, -|- I, and 1^-1-2. The nonbonded inter- 
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depends on the distance between the residues, and 
on their type, cr^ ^ — A,B. The long-range behavior is at- 
tractive for like pairs of residues [C{A, A) = 1, C{B, B) = 
0.5] and repulsive else [C{A, B) = C{B, A) = -0.5]. The 
lengths of all virtual peptide bonds are set to unity. 

In our aggregation study, we have performed mul- 
ticanonical simulations [l5| for two identical pep- 
tides with 13 monomers and the sequence 13.1: 
AB2AB2ABAB2AB is arbitrarily chosen from the Fi- 
bonacci series [13]. For consistency, the simulations were 
repeated for pairs of identical homopolymers, 2 x A13 
and 2 x B13, as well as the larger aggregates 3 x 13.1 and 
4x13.1 [1^ . In all cases, aggregation behaviors of similar 
type as for the two-peptide system were identified. 

For the simulations of the 2 x 13.1 system, the pep- 
tides were confined in a periodic cube with edge lengths 
L = 40. We varied the edge lengths to make sure that 
effects due to this confinement are negligible. A sequence 
of spherical-cap updates 17[ and three-monomer corner 
rotations ensured an ergodic scan of the conformational 
space. After performing 180 multicanonical recursions, 
a total number of 2 x 10^° updates was generated. The 
primary result of these simulations is, up to an unimpor- 
tant constant, the density of states g{E) which has been 
precisely estimated over about 100 orders of magnitude. 

In Fig. [DJa), the microcanonical entropy S{E) — 
In g{E) [ks = 1) flS*] is plotted (up to an unimportant ad- 
ditive constant) for the two-peptide system, ranging from 
the aggregate phase including the lowest energy found in 

the simulation {E^i^ = £^AB,,nin + -E^AB^min + -E^iaLin ~ 

— 18.407), to the phase of the fragmented polymers. The 
conformation of the lowest-energy aggregate has a two- 
cap-like, globular shape with a compact hydrophobic core 
jointly formed by the two heteropolymers, see the in- 
set of Fig. [IJa). It should be noted that the individual 
conformations in the aggregate strongly differ from the 
single-peptide ground states {E^^^^ ~ —4.967 [13]) and 
their respective energies in the aggregate are much larger 
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-3.798). The strongest 



contribution is due to the interaction between the het- 
eropolymers (4Blin ~ -11.412). ^ 

The most interesting region in Fig. [ija) is the phase 



coexistence regime E; 
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-8.85 < E < 1.05 « Efr: 



where the entropy exhibits a convex intruder. The con- 
cave huU Hs{E) = S{Ei,gg) + £^/Tagg, which is the tan- 
gent connecting S'(i?agg) and ^{-Efrag), is the Gibbs con- 
struction. Its slope defines the inverse of the aggrega- 
tion temperature Tagg « 0.198. The interval AQ — 

Sfrag - Sagg = Tagg[^(Sfrag) - ^(^agg)] « 9.90 is the 

latent heat required to release inter-chain contacts at the 
aggregation temperature Tagg. The energy, where the 
difference AS{E) = Hs^E) - S{E) is maximal, is de- 
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FIG. 1: Aggregation transition from the microcanonical per- 
spective: (a) microcanonical entropy S{E) (up to a con- 
stant) and concave hull Tls{E), (b) inverse caloric temper- 
ature T~^{E), and (c) specific heat Cv{E). The errors are 
very small and therefore only shown for Cv{E). 



noted as Esep and the associated maximum deviation is 
the surface entropy A^surf = AS{Esep)- The deriva- 
tive of the Gibbs construction gives the Maxwell line 
Tg^g = const ~ 5.043 in the reciprocal caloric T~^{E) 
curve which is shown in Fig. [Hb). A bijective mapping 
between T and E is only possible for T > T> k. 0.231 
and T < 0.169. This means, for values above T> 

and below r<, that the temperature T is a useful con- 
trol parameter. The two-heteropolymer system forms an 
aggregate for T < T<, where the separation into indi- 
vidual polymers is not useful because inter-polymer at- 
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traction dominates over intrinsic structure formation and 
the aggregate determines the mesoscopic length and en- 
ergy scale. On the other hand, for T > Ty the polymers 
are only weakly interacting fragments, i.e., they can be 
considered separately, the total system energy is an ex- 
tensive variable [E ^-ab ^ab)- Only in these two 
temperature regions, the interpretation of the canonical 
formalism is generic. 

In the transition region T< < T < T>, however, the 
interaction strength between the polymers is as strong 
as intrinsic monomer-monomer attraction and cannot be 
neglected. As a consequence of the convexity of S{E) 
in the interval i?agg < E < i?frag, there is no one-to-one 
correspondence between temperature and energy in the 
transition regime which results in the backbending ef- 
fect: Fragmentation of the aggregate leads to a decrease 
of temperature, although the system energy increases. 
The areas A+ = T~^liE,,p - Ef,,^) - [SiEi,,g) - S{E,,p)] 
and A. = r-i(£;scp-Sagg)-[5(£;scp)-5(i^agg)] formed 
by the MaxweU line and the T~^{E) curve as shown in 
Fig.[ljb) are identical. These areas determine the interfa- 
cial entropy AS'gurt — A^ = A_ which is interpreted 
as the loss of entropy due to the existence of the phase 
boundary [l^l between the aggregate and the fragment 
macrostates of the polymers. Consequently, as the en- 
ergy of the total system is not extensive in the transition 
region, E is the favored control parameter compared with 
T. Therefore, the aggregation transition is more favor- 
ably analysed in the microcanonical ensemble, at least for 
such finite systems like the heteropolymers in our study, 
where an extension towards the thermodynamic limit is 
not possible. 

The most remarkable result is the negativity of the 
specific heat of the system in the backbending region, as 
shown in Fig. [TJc). A negative specific heat in the phase 
separation regime is due to the nonextensitivity of the 
energy of the two subsystems resulting from the inter- 
action between the polymers. "Heating" a large aggre- 
gate would lead to the stretching of monomer- monomer 
contact distances, i.e., the potential energy of an exem- 
plified pair of monomers increases, while kinetic energy 
and, therefore, temperature remain widely constant. In 
a comparatively small aggregate, additional energy leads 
to cooperative rearrangements of monomers in the aggre- 
gate in order to reduce surface tension, i.e, the formation 
of molten globular aggregates is suppressed. In conse- 
quence, kinetic energy is transfered into potential energy 
and the temperature decreases. In this regime, the aggre- 
gate becomes colder, although the total energy increases. 

Figure ^a.) shows the typical bimodal canonical en- 
ergy distribution H{E) ^ g{E) exp(— S/fcsTagg) close to 
the transition temperature Tagg. The maximum points 
are identical with the energies of the phase boundaries, 
i?agg and iSfrag , and the minimum is found at Egep • For 
this reason, the difference of the energies belonging to the 
maximum points of the canoncial distribution is identical 
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FIG. 2: Bimodal canonical energy distribution close to the (a) 
aggregation temperature Tagg and (b) the subphase transition 
near Tagg, 2. Vertical dashed lines mark extremal points. The 
corresponding energies in (a) are identical with those indicat- 
ing the phase boundaries in Fig. [T] and, in particular, to the 
crossing points of the reciprocal caloric curve T~^{E) with 
the Maxwell line in Fig. [TJb). 

with the latent heat AQ. The minimum of this distribu- 
tion coincides with the energy i?scp ~ —3.15, where the 
Maxwell line crosses the T^^{E) curve in the backbend- 
ing regime in Fig. \V[h) . These identifications are easily 
proven by setting the logarithmic derivative of H{E) at 
Tagg to zero, which yields dS{E)/dE = T~^^. The left- 
hand side is the reciprocal microcanonical temperature 
and thus T-^{E) = T^g^. As is seen from Fig. Wih), 
this equation has three solutions, at £'agg, £'sop, and 
-Efrag- Therefore, H{E) possesses three extremal points 
at exactly these energies. Another expected result of 
the correspondence between the canonical and microcani- 
cal representations is that the interfacial surface entropy 
can be written as: AS'gurf = ln(i?(£^agg)/-ff (-^scp)) = 
/cb (£'fiag)/^^(£^sop)) 0]- These expressions serve as 
convenient estimators of the surface tension, which can 
be defined as cr = TaggAS'surfZ-Ragg: where i?agg is the 
square radius of gyration of the aggregate. 

A short remark shall also be devoted to a second, much 
weaker transition that accompanies the aggregation tran- 
sition. It is also of "backbending" type and can be ob- 
served in the fragmentation region in Figs. [Hb) and (c) 
close to i? « —0.32. The associated transition tempera- 
ture is T'agg_2 ~ 0.178 and is, therefore, smaller than Tagg, 
but happens in the energetic region, where the popula- 



4 



tion of fragmented macrostates dominates. In fact, this 
effect is difficult to understand and requires a system 
parameter that allows the structural discrimination be- 
tween macrostates. A detailed microcanonical analysis of 
the square relative distance between the centers of masses 
of the polymers reveals [31 that for energies close to i^frag 
the system is in a fragmented state, and the population 
of aggregated polymers in this energy region is extremely 
small. The situation is different for energies E < 0.22, 
where weakly stable aggregated conformations and poly- 
mer fragments coexist. Only for much smaller energies 
{E < -Eagg), compact aggregates dominate. Having this 
in mind, the transition can also be understood from the 
canonical view. For temperatures below Tagg.2 ~ 0.178, 
stable aggregates (solids) of low energies {E < E^gg) 
dominate. Approaching Tagg,2, the system enters the 
subphase of coexisting unstable pre-molten aggregates of 
comparatively high energies (E « —0.32) and already 
fragmented peptides. From Fig. [TJb) we see that this 
process is also accompanied by cooling due to monomer 
arrangements reducing surface tension. These monomer 
translocations are, however, energetically unfavorable. 
Eventually, for temperatures larger than Tagg, confor- 
mations of weakly coupled separate fragments (liquid) 
dominate. The intermediary subphase is never dominat- 
ing, and therefore unstable. After these remarks this 
transition is already visible in the cusp-like behavior of 
the left, inner wing of the right fragmentation peak in 
Fig. [IKa) . Reweighting to the subphase transition tem- 
perature Tagg.2, the bimodal structure of the canonical 
energy distribution in this energy range is clearly revealed 
in Fig. [2Kb). Compared with the distribution at the ag- 
gregation transition in Fig. [2)Ja), the ratio between max- 
imum and minimum is small and, therefore, also the sur- 
face tension. In consequence, the transition between the 
solid and the pre-molten, unstable aggregates is, com- 
pared with the aggregation transition, negligibly weak. 
Note that the aggregation peak, not shown in Fig.l^Jb), is 
much more pronounced than the peaks of the pre-molten 
aggregates at E w —0.32 and fragments a.t E 0.73. 

In this Letter, we have shown by employing a 
mesoscopic hydrophobic-polar heteropolymer aggrega- 
tion model that the aggregation transition is a phase sep- 
aration process, where the loss of entropy due to the ex- 
istence of the phase boundary results in negative specific 
heat. This is an effect which is guided by changes of the 
interfacial entropy as a result of surface effects. There- 
fore, this effect is expected to disappear in the thermody- 
namic limit of macroscopic systems. It should strongly 
be emphasized, however, that peptides and proteins, like 
the exemplified model heteropolymers used in our study, 
are necessarily systems of finite length and a thermody- 
namic limit cannot be defined. For this reason, standard 
canonical formalisms for the analysis of conformational 
pseudophase transitions with phase separation are not 
suitable for these systems, since the temperature is not 



a unique control parameter and the total system energy 
measured in units of energy scales of mesoscopic particles 
(e.g., aggregates or single polymers) is not an extensive, 
separable quantity. In such cases, microcanonical ther- 
modynamics with the energy itself as the external con- 
trol parameter provides a more favorable basis for the 
study of first-order- like transitions. The interesting phe- 
nomenon of the negativity of the microcanonical specific 
heat in peptide aggregation should be motivation for an 
experimental verification which is still pending. 
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